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Residual stress. 
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One-fifth of tensile yield stress. 
Tensile yield stress. 

Test identifying information. 
Variable equivalent to SIGMIN. 


Variable equivalent to SIGMAX. 


12 





UN TRODUCTION 


In an aircraft fatigue life monitoring program, deter- 
mination of the actual loading environment which produces 
the fatigue is a major oplem. With both military trainer 
and fighter aircraft, which may be used for a variety of 
duties, regular or continuous load recording programs have 
to be considered mandatory. Such load recording programs 
provide the data to calculate the consumed fatigue life of 
Operating aircraft and are used in selecting design spectra 
for future aircraft designs. 

flor facilitate the calculations of fatigue life from 
such data, it is important to know not only the loading 
magnitudes but also the time sequence of the loading. The 
igue life of a local critical point varies due to 
residual stresses remaining after the application of a load 
@ausing local plasticity in tension or compression. For 
instance, a peak tensile load into the plastic flow region 
leaves a residual compressive stress, which lowers the 
local magnitude of the following tensile stresses, thus 
increasing the fatigue life. In a similar manner a peak 
compressive load into the plastic flow region, if it were 
to occur, would leave a residual tensile stress, which 
decreases the fatigue life. 

Currently, fatigue monitoring of naval aircraft is 


based on the total number of g readings recorded at four 


ES 





selected levels by an exceedence level counting accelerometer. 
Using micropressors in the near future it will be possible 

to record each maximum and minimum load level experienced 

by an aircraft in sequence. The data collected can be used 

to monitor the fatigue life of a structure via the deter- 
mination of damage accumulated at a point found to be 
Simre1cal in a structural test of a prototype. It should 

be noted that this calculation uses a theoretical model. 
Damage is not observable or measurable. 

The objective of this thesis is to use a computer program, 
employing the range-pair counting technique on time sequence 
recorded maximum and minimum aircraft loads and a relaxation 
usdel to consider residual stresses, to calculate plastic, 


estic, and total aircraft structural damage. 
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ПОЕ Е ЕО DESCRIPTION. OF THE COMPUTER PROGRAM 


The computer program used in this fatigue analysis is 
divided into four modules. Each module is clearly labeled 
in the program listing (page 93). In module I the input 
data are read and assimilated in preparation for future 
calculations. The input parameter requirements are ου 
in detail beginning on page 18, and six sample sets of 
input data are illustrated starting on page 40. 

When the load sequence randomization technique is used 
input loading data, it is a part of module I. The load 
sequence randomization technique uses the computer library 
subroutine RANDU to place 10 percent of the MIL spectrum 
A positive loads (Table I) in a random order. Each of the 
positive loads is paired with a minimum load of 11 percent 
limit load, or l-g. During the randomization each load has 
am equal probability of selection. <A counter restricts the 
number of times a value is selected to the number of 
occurrences of the particular load level in MIL spectrum A. 

Nova tae tons оп the randomization of the MIL spectrum 
A load levels are available. The first option permits 
randomization of the negative MIL spectrum A load levels 
using the same technique as in the case of the positive 
MIL spectrum A load levels. Each negative load level is 
randomly paired with a positive load level. Since there 


are not as many negative load levels in MIL spectrum A as 
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Tree I 


FREQUENCY OF MANEUVER LOADS 


NUMBER OF TIMES PER THOUSAND HOURS THAT 
TODA ACTOR Т5 EXPERIENCED 


PERCENT OF MAXIMUM (POSITIVE) FLIGHT MANEUVER LOAD 
SYMMETRICAL LIMIT LOAD FACTOR SPECTRUM A 
55 17000 
45 9500 
35 6500 
65 4500 
ТЭ 2500 
85 1360 
3 440 
105 150 
115 40 
175 16 
TOTAL 42006 
PERCENT OF MAXIMUM (NEGATIVE) FLIGHT MANEUVER LOAD 
EREMETRICAL LIMIT LOAD FACTOR SPECTRUM A 
0 500 
10 200 
20 100 
50. 60 
40 33 
50 А 50 
60 25 
70 20 
80 15 
90 10 
100 5 
110 3 


TOTAL 1003 
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there are positive load levels, the excess load levels are 
paired with l-g loads. Mixing of the negative MIL spectrum 
A load levels and the l-g load levels is accomplished using 
the computer library subroutine RANDU. For this ordering 
thereis an 80 percent chance that each negative load selected 
will be a l-g load. A counter ensures the proper numbers of 
each load type are included. Eighty percent probability of 
selection of a l-g load was used to spread the smaller 
number of negative MIL spectrum A loads throughout the 
sequence, and also to speed computer operation by not having 
the matrix of negative MIL spectrum A loads addressed so 
often, after the counter indicated the matrix elements had 
all been used. 

The second option for randomizing loads provides for the 
inclusion of a ground cycle between each flight. The ground 
cycle loading is taken as -2500 psi, or -8 percent limit 
load. Eight percent is the level of the ground cycle load 
used by the A7 manufacturer. There are 4201 positive 
loading events in 100 hours of MIL spectrum A flight time. 
Ber flieht is arbitrarily taken as one hour in length, so 
a ground load is paired with every forty-second positive 
loading event. 

In module II the local stresses and strains are deter- 
mined from the input data. If the local peak stress is in 
the plastic region, the specimen is assumed to unload 
Ei tccally, leaving a local residual compressive stress 


in the material. Іп the analysis developed by Potter [Ref. 6], 
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the transient portion of the residual stress so produced 
relaxes toward zero, or an equilibrium residual stress. 
After the stresses and strains are calculated, the local 
stress cycles are counted in module III using the range- 
pair counting technique. In module IV the fatigue damage 
EN !culated. Damage 1s determined separately for elastic 
and plastic strain events. The elastic and plastic damage 
tor each flight is listed as output, as well as the total 


cumulative damage from any previous flights. 


A. INPUT DATA REQUIREMENTS (MODULE I) 

Samples of data are listed starting on page ho; format- 
ting is described beginning on page 93. Both pages should 
be consulted while reading the following descriptions. 

E Data Card 1 

Only one of these cards is required or permitted 
for each program run. 

NDECK = the number of data decks to be run sequen- 
ШОШ у The input on this card will be common to all data 
decks run. It is not necessary for the different data 
decks to have any parameters in common. 

ORNIS ἵπμο να πο controlling the write statements. 
Output available during the process of the analysis 
includes: 

a. Maximum and minimum applied stress of each cycle 
and the local stress response throughout the spectrum. Also 
printed out is the residual stress, equilibrium stress, 


ПЕРИ са cycles, and the equilibrium period for each cycle. 
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ο ο τ]ς local stress history as input into 
the range-pair subroutine and the resulting range-paired 
spectrum. 

сах plastic Strain occurrence during the 
Spectrum and the damage associated with each strain reversal. 

d. The accumulated damage associated with the plastic 
strains. 

e. The range-paired elastic stress spectrum and the 
damage associated with each level for block loading, or for 
each cycle for RANDU generated loadings. 

f. The accumulated damage associated with the cur- 
rent block of loading, including the plastic strain damage, 
and the total damage since the initiation of cycling. 

πμ ше ass tened to ITPRINT controls which pieces 
data are included in the printout. 

If IPRINT = 1, all six items listed above are 
printed for each flight or block of loads. 


ΠΡ 2, all items listed above except b, 


are printed. 


If IPRINT 3, only items d and f listed above are 
printed. : 

An example of each output option is given beginning 
on page . When the number of loads is more than fifty, 
the output is auite voluminous, and a paper usage default 


may interrupt program operation. With large numbers of 


loads, IPRINT should be set equal to 5. 





ПЕРЕМ = the value controlling entry into the range- 
paar counting subroutine. 


AR] CM 


Кеше nhange-palr Counting méthod 15 used. 


Tí DOM 


пе \тапое-ратг counting method is 
skipped. 

Each data deck run sequentially must contain the 
following cards. The number of data decks was indicated by 
the value assigned NDECK. 

ee Card 1 

Information on this card is descriptive, alphanumeric 

data used only as a heading test identification. 
p card 2 

The name of the material type and four material con- 
stants are listed on this card. The material type is an 
alphanumeric entry. The tensile yield stress, LCF strain 
intercept, inverse of the Coffin-Manson slope, and the 
modulus of elasticity are real number entries. The LCF 
strain intercept is the numerical value of the ordinate 
intercept in the log-log Coffin-Manson plot of plastic 
Tin range vs. cycles to failure. Units of tensile yield 
stress and modulus of elasticity are ksi. 

Cards 5 Through 6 

Each card contains the three coefficients of a second 
order least-squares curve fit for S-N data obtained from a 
Goodman diagram published in MIL-HDBK-5A. The curves are 
Ent for lives of E 10?, ο... and 10” cycles. Space is 


provided on each card for an alphanumeric indication of 
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'. πια source. Only the data source printed on card 6 
LU be retained in Computer storage for subsequent printout. 
uns 72 through 80 are used for information only and are 
not read by the computer. 
oe Card 7 
Constants to be used in the equilibrium period cal- 
culations are listed on this card. 
oe Card 8 
Dhcgscressmeoncentration factor is the only entry 
on this card. 
fee Card 9 
The first entry on this card is the number of blocks 
of loads, or the number of flights. The terms block and 
flight are used interchangeably. This entry is the number 
of times the list of loads is to be repeated. The next 
entry is the number of loads. The number of loads will 
always be the same as the quantity of cards 11 used as input. 
The last entry is the number of types of loads. Loads are 
indicated as being a different type in order to alter the 
loading pattern during subsequent flights. It is important 
to understand that the program will always consider load 
types in ascending numerical order. A flight listing load 
types in the order 1, 2, 3, 4 will produce identical results 
КЕ а рне listing load types in the order, 4, 3, 2, l. 
8. Card 10 
Only the limit load is listed here, with units of 
ο. The maxima and minima of all cycles are input as a 
decimal fraction of this load entry. 


aL 





o Cards 11 

Information identifying each load cycle is presented 
on these cards. The first entry is an integer representing 
the step, or number of the load. This integer increments by 
one on each card, with the last card having a value equal to 
ШТЕУЕІ,. Other entries on these cards are the type of load, 
the minimum and maximum decimal fraction of the limit load 
for the particular step, and the number of consecutive 
cycles at this load level. It is important to notice that 
the format requires the number of cycles to be entered as 
a real number, although fractions of cycles will not be 
processed. 

ШЕ Cards NLEVEL + 11 

The first entry 1S an integer representing the number 
of the block, or flight. This entry is a dummy variable and 
is not used in any subsequent calculations. The other 
entries on this card indicate which types of loads are to 
be processed on this particular flight. Each type of load 
may be processed only once on each flight. The order in 
which the types are listed is of no significance, since the 
load types are automatically processed in ascending numerical 


order. 


B. RANGE-PAIR COUNTING METHOD 

In assessing the life consumed by an individual aircraft, 
damage is calculated from Miner's Law using S-N data. This 
ШІӘЕ of life calculation is of greater value if the cycle 


counting method used takes into account the actual load-time 
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history of the aircraft. Many load cycle counting methods 
have been developed and used. Fatigue test experience indi- 
ατα useful counting technique for aircraft fatigue 
must take into account the loading sequence as well as the 
magnitude of the maximum and minimum load peaks [Ref. 10]. 
The two most realistic counting techniques are the range- 
pair counting method and the rain-flow counting method. 

Both methods experimentally yield approximately the same 
results. The major difference between the two methods is 
that the rain-flow method counts in terms of load ranges, 

or half cycles, whereas the range-pair method records the 
time history in terms of complete load cycles. When counting 
in terms of full cycles, the two methods are equivalent. 

The merit of counting half cycles is not important when 
analyzing aircraft fatigue because the large number of load 
reversals does not permit a half cycle to influence the 
numerical results. The range-pair counting method corresponds 
to the stable cyclic stress-strain behavior of a material 

in that strain ranges counted as cycles will form closed 

E -strain hysterisis loops [Ref. 4]. 

The computer program employed in this thesis has the 
option of using the range-pair counting method, or simply 
using each peak in the order that it occurs. When the range- 
pair counting method is not used, only momentary load values 
are considered, and information regarding the cyclic stress- 
rin pattern, which is important in fatigue calculations, 


EN t Also, minor load variations, which"are of little 
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шівоттапсе іп fatigue calculations, are counted as additional 
e cles. 

Mer anee pair counting method is illustrated in Figure 1. 
It counts a strain range as a cycle if it can be paired with 
a subsequent straining of equal magnitude in the opposite 
direction. For a complicated load history, some of the 
ranges counted as cycles will be simple ranges, such as 2-3, 
during which the strain does not change direction, but others, 
such as 1-8, will be interrupted by smaller ranges which will 
also be counted as cycles. In Figure 1 ranges are marked with 
solid lines and the paired ranges with dashed lines. 

Each peak is taken in order as the initial peak of a 
range, except that a peak is skipped if the part of the 
history immediately following it has already been paired with 
a previously counted range. If the initial peak of a range 
is a minimum, a cycle is counted between this minimum and 
the most positive maximum which occurs before the strain 
becomes more negative than the initial peak of the range. 

For example, in Figure 1 a cycle is counted between peak 1 
and peak 8, peak 8 being the most positive maximum that 
occurs before the strain becomes more negative than peak l. 
If the initial peak of a range is a maximum, a cycle is 
counted between this maximum and the most negative minimum 
which occurs before the strain becomes more positive than 
the initial peak of the range. For example, in Figure l a 
cycle is counted between peak 2 and peak 3, peak 3 being 


the most negative minimum before the strain becomes more 
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Верите 1 


Illustration of Range-Pair Counting Method 
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positive than peak 2. Each range that is counted is paired 
with the next straining of equal magnitude in the opposite 
direction, explaining why complete cycle rather than half 
cycle counts are made. For example, in Figure 1 part of the 
range between peaks 8 and 9 is paired with the range counted 


between peaks 1 and 8. 


C. DAMAGE CALCULATION 

Damage due to plastic strain and damage due to elastic 
strain is calculated separately. The behavior of the specimen 
memassumed to be elastic-perfectly plastic as illustrated in 
Figure 2. Behavior in compression is assumed to be the same 
as the behavior in tension. The total damage is determined 


using the Miner's summation D = X£Xn/N where 


D total damage 
n = number of cycles at a particular load level 


N = number of cycles to failure at a particular 
load level 


Failure of the specimen is defined as occuring when D equals 
pnity. 

Etic strarzn damage is calculated using the Coffin- 
Manson theory [Ref. 5]. The Coffin-Manson theory is based 
on an experimentally determined log-log plot of plastic 
Strain range vs. cycles to failure as shown in Figure 5. 

The slope of the line in Figure 3 is approximately  -1/2 
for all metals, with the ordinate intercept making the 
behavior of each metal unique. The equation used for cal- 


Sulation of plastic strain damage is 


26 





stress, 6 


tensile yield EN ” 
stress, 6 ty 


strain, E 


Figure 2. ELASTIC-PERFECTLY PLASTIC 
BEHAVIOR 






moe plastic | e 
strain range, Slope approximately -3 


ED 


log cycles to failure, N 


Figure 3. COFFIN-MANSON DIAGRAM 
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c NAE 


where N number of cycles to failure at a specific plastic 


strain 

e, ΙΟ: the plastic strain 

k = slope of the Coffin-Manson curve, approximately 
- 1/2 


c = ordinate intercept; experimentally for each metal 
type 

For strain levels in the plastic range, failure occurs in 
less than ο. cycles, and damage is calculated using the 
Coffin-Manson theory. For strain levels in the elastic range, 
the Goodman diagram is the basis for calculation of damage. 
The maximum and minimum local strain levels are sequentially 
compared with second order least-squares curve fit data from 
the Goodman diagram. The curve fits used are for lives of 


7. 10°, and 10’ 


cycles [Ref. 3]. The elastic damage 
is totaled using the Miner summation and combined with the 
meeastic damage, computed in a similar way, to yield the total 
damage. The amount of plastic damage and the total damage 


for a particular flight, and the total damage including all 


prevrous flights is printed at the end of each flight. 
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ο CUSSION ОЕ RESULTS 


The load spectrum used in this analysis is the load spectrum 
A from MIL-HDBK-5A, and the loads are assumed to be applied at 
nos tation 32 on an A7 aircraft wing [Ref. 2]. Load data 
were input as block loads and also as single randomized cycles. 
The Miner's Law summation of fatigue damage for this loading 
spectrum is .1300 per 1000 flight hours. Results obtained 
using this computer program are given in Tables II, III, and 
mee it is important to notice that when fatigue damage is 
GC ulated without using the range-pair counting method, it 
makes only a minor difference on the total fatigue damage when 
the sequence of loading is changed or when negative loads are 
included. 

ΙΕΚ loads were arranged in LO-HI, HI-LO, and HI-LO-HI 
E:uences. The results indicate the LO-HI sequence to be 
more damaging than the HI-LO sequence, with damage due to 
the HI-LO-HI sequence falling between the other two block 
load sequences. In the HI-LO sequence the HI loads leave 
a local residual compressive stress in the material, and 
the LO loads do less damage. 

In the case of RANDU generated load sequences, identical 
starter integers were used in order to evaluate the effect 
on the damage of inclusion of MIL spectrum A negative loads 
and ground cycles. The more negative loads that are included, 
the worse is the fatigue damage. The increase in damage with 


more negative loads is caused by the negative loads decreasing 
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TABLE II 


FATIGUE DAMAGE DUE TO POSITIVE MILSPEC A 
BLOCK LOADING, NORMALIZED TO 1000 HOURS 


SEQUENCE OF 


BLOCK LOADS ΕΡΕ NOT USINGIRP CMe 
ИШЕН! aro 30102 „16780108 
> LO . 069245014 .069243014 
ES LO - A 1 ‚10864216 ος ENS 


* RANGE-PAIR COUNTING METHOD 


TABLE III 


FATIGUE DAMAGE DUE TO MILSPEC A 


RANDOMIZED LOAD INPUTS, USING RPCM, NORMALIZED TO 1000 HOURS 


EUEGERS USED 


TO START RANDU MINIMUM LOAD VALUES 

iene EN i NEGATIVE GROUND 

LIMIT LOAD MICSPEC A CX GUES 
83745,54711,54487 05922726558 “020267735 .060848035 
15547,66549,7 .095124505 03530651607 10429338 
 5955,4777 “0455321515 .043640956 .036781987 
48621,3,491 .018672124 054077277 029881012 
73,559,1001 0513919174 .042006299 2057120679 
E 855,1 .037665060 .049336702 050375555 
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TABLE IV 


FATIGUE DAMAGE DUE TO MILSPEC A RANDOMIZED LOAD 


INPUTS, NOT USING RPCM, NORMALIZED TO 1000 HOURS 


ШҮРЕСЕК5 USED 
TO START RANDU 


iis PERCENT 

LIMIT LOAD 
83745,54711,54487 .040602274 
547, 0060549,7 . 040157805 
0556835, 4777 .041204450 
48621,3,491 ПОЧ 60218 
55559,1001 .040792227 
B5 335,1 .041056550 


31 


MINIMUM LOAD VALUES 


NEGATIVE 
MILSPEC A 


20375393917 
047727041 
.048346408 
1049250121 
20470453538 


.049220584 


GROUND 
CYCLES 


.049124956 
.049612150 
.049540661 
-050957575 
04.931.966 
.050481893 





the local residual compressive stress caused by positive loads 
miro tne plastic zone. This permits the local positve elastic 
stresses to be more damaging. 

Six plots are presented, beginning on page 33, to illus- 
Exe the pattern of the individual load randomization process. 
ШЕШ рІот рогтгтаув the first two flight hours of the aircraft. 
The range of values in the columns of Table III are caused by 
such variations in the load randomization pattern. 

The fatigue damage resulting from randomization of loads 
on a cycle basis, as shown in Table III, is significantly 
 еЕтетеп: from the damage due to the HI-LO-HI block loads 
illustrated in Table II. Vought Aeronautics Division con- 

n cd a fatigue study using randomized flights. Each flight 
was based on a mission profile specified by the Air Force. 
Fatigue damage resulting from the randomized flights was in 
Grose agreement with fatigue damage resulting from HI-LO-HI 


block loading. 
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IV. RECOMMENDATIONS FOR FUTURE RESEARCH 


Subsequent work in this area should be directed toward 
refining the input load sequence. The Naval Air Development 
Center collects accelerometer data from individual aircraft 
on a monthly basis. These load data could be used in this 
computer program to see if randomizing the load sequence over 
the life of the structure provides a different picture of the 
fatigue damage than randomizing by the month. Also, research 
is needed to determine the magnitude and sequence of the 
ground cycles. In this analysis all ground cycles were con- 
Sidered of equal magnitude and equal spacing, which is not 
the case with a real airplane. A landing load spectrum should 
be used to be more accurate. 

The computer program itself could be improved by a better 
understanding and definition of the relaxation phenomena. 
Experiments to quantify the relaxation parameters will enhance 
the program. 

The Air Force is monitoring crack growth as the basis for 
their fatigue damage assessment. It would be productive to 
compare the total fatigue damage calculated by the two 


methods. 
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